Statistical properties of the quantum anharmonic oscillator in one spatial dimension. 
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1 Abstract 

The random matrix ensembles (RME) of Hamiltonian matrices, e.g. Gaussian random matrix 
ensembles (GRME) and Ginibre random matrix ensembles (Ginibre RME), are applicable 
to following quantum statistical systems: nuclear systems, molecular systems, condensed 
phase systems, disordered systems, and two-dimensional electron systems (Wigner-Dyson 
electrostatic analogy). A family of quantum anharmonic oscillators in one spatial dimension 
is studied and the numerical investigation of their eigenenergies is presented. The statisti- 
cal properties of the calculated eigenenergies are compared with the theoretical predictions 
inferred from the random matrix theory. Conclusions are derived. 

2 Quantum harmonic oscillator in. D = 1 spatial dimen- 
sion. 

Firstly: Let us consider the Hilbert space: 

Vi = L2(R,C,dx), (1) 

of the complex valued wave functions that are (modulus) square integrable on the set 
R of the real numbers, and C is the set of the complex numbers. The Hilbert space Vi is 
separable space, and its orthonormal basis B is a set of Hermite's functions \E'r (Fok's func- 
tions, eigenf unctions of the dimensionless Hamiltonian operator 7i of the quantum harmonic 
oscillator in D = 1 spatial dimension): 

^,(a;) = NrHr{x)exp{-^x^),Nr = [v^r!2^]5,r e N, (2) 
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where N is a set of natural numbers including zero, whereas: 

Hr{x) = {-lYeM^')-^^M-^'), (3) 

is rth Hermite's polynomial. We assume from now that all the operators are dimensionless 
(nondimensional) . The dimensionless quantum operator x of the xth component of the 
position vector (radius vector) acts on the Hermite's basis function as follows: 



X^r = + y (4) 

whereas the quantum operator of the xth component of linear momentum vector is given 
in the basis by the following formula: 



T 1 f^T 1 /r + 1 

Px*r= T^-^'r-l-Ty^-^r+l- (5) 

Matrix elements x^r and {Px)i,r oi these operators equal correspondingly: 



pF r -\- 1 

Xl,r = (^il^^r) = \ 7:\r-l + \ —rS^r+l, (6) 



and 



where 



1 fV 1 r + 1 

{Px)l,T = (*i|Px^r) = -\l2kr-l - -y 5i,r+l, (7) 

kr = I J; j ^ ^ , (8) 

is discrete Kronecker's delta (it is not continuous Dirac's delta). The quantum operator x^ 
of the square of the xth component of radius vector after acting on the basis function reads: 

x^'^r = \ Wr - lv^*^-2 + (2r + 1) + VrTTv^r + 2*^+2] , (9) 
and its matrix element {x'^)i,r is equal: 



[X' 



\r = -[Vr - lV^-6i,r-2 + (2r + l)6i,r + VV+lVV+26i,r+2\, (10) 



whereas the quantum operator of the square of the xth component of the linear momentum 
vector has the following representation in the basis: 

Pl'^r = h-Vr + (2r + 1)*^ - Vr + Wr + 2*^+2], (11) 



2' 
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and its matrix element {pl)i,r is equal: 

Therefore, the dimensionless (nondimensional) quantum Hamiltonian operator reads: 

n^pl + x', (13) 

acting on the basis function, it produces: 

n^r = (2r + 1)*„ (14) 

and its matrix element 7ii,r is equal: 

ni,r^{2r + l)Si,r, (15) 

where 

= 2r + l, (16) 

is the rth eigenenergy of Ti. The eigenenergies are simply all odd natural numbers, and the 
quantum Hamiltonian is diagonal operator, and its matrix representation is diagonal oo x oo 
matrix. Note, that if one introduces the notion of nearest neighbour energy spacing (NNS): 

Sr = er+i-er,r = l,---{N-l),N>2, (17) 

then for the quantum harmonic oscillator it holds: 

Sr = 2 = const, (18) 

so the eigenenergies are equidistant. Here N is the number of energies dealt with (it is 
also a dimension of truncated Hilbert space Vi,7v = C^). The eigenfunctions of Vi.tv are 
A^-component constant complex vectors (analogs of spinors), and the operators acting on it 
are N x N deterministic complex-valued matrices. The probability distribution P/v-i of the 
spacings is discrete one point distribution for any finite value oi N, N > 2: 

tending in the thermodynamical limit — > oo to the singular Dirac's delta distribution: 

Poo{s)^5{s-2). (20) 

Secondly, let us perform more difficult task consisting of calculating all the absolute 
moments: _ 

/■oo 

{mn)l,r ^rUn^ (x"),,, = {^l\x^'^r) = / n(x)x"*,(x)dx, (21) 
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of the position operator x in the Hermite's (Fok's) basis B. One can calculate the lower 

moments manually, e. g., using recurrence relations, matrix algebra, etc., but it is tedious 
(even for 3 < n < 6). If one wants to calculate all the moments then he must return to the 
beautiful XIX century mathematics methods and after some reasoning he obtains the exact 
formula: 



{mn)l,r ^rUn^ {.x'')l,r = 

^ ^ ^ ^hto ^ M-2j)\k\{r-2ky. 

where [•] is entier (step) function, F is Euler's gamma function. Therefore, the matrix 
representation of the even power operators x'^^ in the basis B are hermitean (symmetrical real) 
matrices with nonzero diagonal and nonzero p subdiagonals (and nonzero p superdiagonals), 
where the distance of the nearest superdiagonals (or subdiagonals) is 2 (the diagonal is also 
distant by 2 from the nearest super- and sub-diagonal), whereas the odd power operators 
^2p+i -j^ ^YiQ basis B are hermitean (symmetrical real) matrices with zero diagonal and p 
nonzero subdiagonals (and p nonzero superdiagonals), where the distance of the nearest 
superdiagonals (or subdiagonals) is 2 (the nearest super- and sub-diagonal are also distant 
by 2). The physical interpretation of the superdiagonals (and subdiagonals) is connected 
with the absorption (emission) of phonons. 



3 Quantum anharmonic oscillator in D = 1 spatial di- 
mension. 

Thirdly, we are ready to deal with the quantum anharmonic oscillator in D = 1 spatial 
dimension. Its dimensionless Hamiltonian operator T^f anharm reads: 

s 

^l,anharm — ^ + dgX^ , (23) 
s=0 

where S* is a degree of the anharmonicity of the oscillator, and the prefactors are the 
strengths of anharmonicity. The matrix elements of the anharmonic Hamiltonian operator 
are: 

s 

('^f,anharm)/,r = ^r^Z.r + Yl «s(£')/,r, (24) 

s=0 

where the representation of the quantum anharmonic oscillator in the quantum harmonic 
oscillator basis B is mathematically correct, because the basis B is a complete set, and the 
Hilbert space of the eigenfunctions of the anharmonic oscillator is isomorphic to the Hilbert 
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space Vi for the harmonic oscillator, provided that the anharmonic potential: 



Wl%a.m(^) = E«.^^ (25) 
s=0 

is bounded from below (there are no scattering eigenstates). 

Fourthly, we repeat the "Bohigas conjecture" that the fluctuations of the spectra of the 
quantum systems that correspond clasically to the chaotic systems generally obey the spectra 
of the Gaussian random matrix ensembles. The quantum integrable systems correspond to 
the classical integrable systems in the semiclassical hmit [1, 2]. The probabihty distributions 
Pff of the nearest neighbour spacing for the Gaussian orthogonal ensemble G0E(2) of 2 x 
2 Gaussian distributed real- valued symmetric random matrix variables {f3 = 1), for the 
Gaussian unitary ensemble GUE(2) of 2 x 2 Gaussian distributed complex- valued hermitean 
random matrix variables (/? = 2), for the Gaussian symplectic ensemble GSE(2) of 2 x 2 
Gaussian distributed quaternion- valued selfdual hermitean random matrix variables {(3 — A), 
and for the Poisson ensemble (PE) of the random diagonal matrices with homogeneously 
distributed eigenvalues on the real axis R are given by the formulae: 

Pp{s) = 9{s)Afssfexp{-Bps^), (26) 

for the Gaussian ensembles, and 

Po(s) = ^(s)exp(-s), (27) 

for the Poisson ensemble, where 9 is Heaviside's unit step function [3, 4, 5, 6, 7, 8, 9, 10, 11, 
12]. The constant numbers: 

_ T^^\{f3 + 2)/2) _ r^((/3 + 2)/2) 

"^'-^f^^W^m) ^^-rw + i)/2) ^ ^ 

are given by the formulae: = 7r/2, Si = 7r/4 (GOE), A2 = 32/7^^ = Vtt (GUE), 
and A4 = 262144/7297r^, B4 = 64/9% (GSE), respectively. For the Gaussian ensembles the 
energies are characterized by the "level repulsion" of degree f3 near the origin (vanishing 
spacing s = 0), and the probability distributions vanish at the origin, and the quantum 
system with "level repulsion" is cast to the class of quantum chaotic systems. For the Poisson 
ensembles the energies are characterized by the "level clustering" near the origin (vanishing 
spacing s = 0), and the probability distributions Pq has maximum at the origin, and the 
quantum system with "level clustering" are treated as the quantum integrable system. After 
many numerical experiments conducted with different quantum anharmonic oscillators (up 
to the sextic quantum anharmonic oscillator 5* = 6) we draw conclusion that majority of 
them behaves like quantum chaotic system, the eigenenergies tend to cluster, the histogram 
of nearest neighbour spacing is closer to the Pq distribution resulting from the Poisson 
ensemble. 
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